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Isentropic thermodynamics in the PNJL model 
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We discuss the isentropic trajectories on the QCD phase diagram in the temperature and the quark 
chemical potential plane using the Nambu-Jona-Lasinio model with the Polyakov loop coupling 
(PNJL model). We impose a constraint on the strange quark chemical potential so that the strange 
quark density is zero, which is the case in the ultra relativistic heavy-ion collisions. We compare 
our numerical results with the truncated estimates by the Taylor expansion in terms of the chemical 
potential to quantify the reliability of the expansion used in the lattice QCD simulation. We finally 
discuss the strange quark chemical potential induced by the strangeness neutrality condition and 
relate it to the ratio of the Polyakov loop and the anti-Polyakov loop. 
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I. INTRODUCTION 

Thermodynamic properties of hot and dense matter 
out of quarks and gluons are of theoretical and experi- 
mental importance. To this end the Monte-Carlo method 
on the lattice has worked quite successfully to simulate 
the matter at high temperature from the first principle, 
that is. Quantum Chromodynamics (QCD) P, [2, H, 0, [1] . 

It is believed that such hot and dense matter has been 
created in the Relativistic Heavy Ion Collider (RHIC) lo- 
cated at BNL. Since the baryon stopping power in the 
nucleus-nucleus collision is small at the RHIC (top) en- 
ergy -^Sjviv — 200 GeV, the created matter is nearly free 
from the net baryon density, and thus, the correspond- 
ing baryon chemical potential is much smaller than the 
temperature. 

We are facing two experimental possibilities for the 
future exploration of the QCD phase diagram: One is 
going toward even hotter matter as planned in the Large 
Hadron CoUider (LHC) at CERN. The other one is re- 
alizing denser matter, which is accessible by collisions at 
smaller energy with larger baryon stopping power. The 
formation of baryon-rich matter is within the scope of the 
Facility for Antiproton and Ion Research (FAIR) at GSI 
and the systematic energy scan is also under discussion 
in the future plan for RHIC with emphasis on the QCD 
critical point search. The present work is focused on 
the latter; the thermodynamic properties of matter with 
a finite quark chemical potential /i whose magnitude is 
comparable to the temperature T. 

The lattice QCD simulation has a serious limitation if 
a finite chemical potential is turned on in the Dirac oper- 
ator @, 0]- That is, the notorious sign problem prevents 
us from applying the Monte-Carlo method to a finite- 
density system. There are a number of proposals to tame 
this problem among which the Taylor expansion in 
terms of ^/T seems to work well insofar as /i/T is un- 
der the radius of convergence [13, [IH ■ In this way the 
equation of state (EoS) at finite T and fi is partially avail- 
able from the direct lattice QCD simulation. Actually the 
EoS is an indispensable input for the sake of the hydro- 
dynamic evolution of matter. From the s/ns-constant 



line, where s is the entropy density and is the baryon 
number density, we can draw the isentropic trajectory 
along which the adiabatic system evolves [lOl, u3L 

An interesting progress has been made in ReL 10[; 
the isentropic thermodynamics was investigated with a 
constraint that the strange quark density must be zero, 
i.e. Us = 0. Such a constraint is necessary to emulate 
matter created by the ultra relativistic heavy-ion colli- 
sions because thermalization is achieved within the time 
scale of the strong interaction and the system is far from 
/3-equilibration. 

In this paper we will utilize the Nambu-Jona- 
Lasinio model with the Polyakov loop coupling (PNJL 
model) [11, M, [H, [13, [S El H mi to reveal the isen- 
tropic thermodynamics with = imposed. In contrast 
to the NJL model study along the same line [l^j (see also 
Ref. Q for another type of approach) , the PNJL model 
has an advantage that a part of the gluon degrees of free- 
dom is included, which gives us a hope that the PNJL 
model can lead to a better EoS refiecting the gluody- 
namics as well as the chiral dynamics. In Ref. |24| we 
can already find the isentropic trajectory evaluated in 
the PNJL model and notice that they are qualitatively 
similar to the results in Ref. [22i|. We will see that the 
constraint Us ~ causes only little change in the resul- 
tant isentropic trajectory because the strange quark mass 
Ms is still substantially large around the crossover region 
and the strange density is suppressed in any case. 

The extension from the NJL model to the PNJL model 
is not a minor improvement, however. Our nontrivial 
finding is, in fact, that the PNJL model is capable of cap- 
turing the correct behavior of the induced strange quark 
chemical potential /i^ > to keep Us = 0. Since the 
nonzero /i^ has its origin in confinement physics, as we 
will discuss later, the NJL model without any confine- 
ment effect is of no use but the PNJL model naturally 
provides us with /i^ > which is related to the ratio of 
the Polyakov loop i and the anti-Polyakov loop £. As a 
matter of fact, although it is known that £ deviates from 
£ at ^ 7^ [§, [13, m, , our present work is the very 
first demonstration to show that this difference £ i has 
a physical consequence in the best of our knowledge. 

This paper is organized as follows: We explain the 
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model definitions in Sec. |TT1 Then, in Sec. IIIII we show 
the numerical results with the constraint Ug — 0, followed 
by the discussions on the validity of the Taylor expansion 
in Sec. IIVI We elucidate the physical meaning of /ig 7^ 
in Sec. Mand the summary is in Sec. |Vll 



II. MODEL SETUP 

The thermodynamic potential of the three-flavor PNJL 
model consists of three pieces, 0. = r^vacuum + ^^thormai + 
f^Poiyakov, namely, the vacuum part (zero-point energy 
and the condensation energy), the thermal part, and the 
Polyakov loop potential, respectively. We shall take a 
close look at them in order. 

The vacuum part is exactly the same as in the ordinary 
NJL model; 



-6E 



' d^p 
(2^ 



gs(^{uu)^+ {ddf+ (ss)^] + Agu{uu) (dd) (ss) , 



(1) 



where i refers to the quark flavor running over m, d, and 
s. The energy dispersion relations are ei{p) = \/ + Mf 
with Mi = mi - 2gs{qiqi) - 2gueijk{qjqj){qkqk), where 
gs represents the four-Fermi coupling constant and 
represents the 't Hooft interaction strength. There are 
three more model parameters: the light current quark 
mass ruu = rud, the heavy current quark mass rus, and 
the ultraviolet cutoff A. We here adopt the parameter 
set by Hatsuda-Kunihiro [27l |: 

A = 631.4MeV, 

gsA' = 3.67, goA^ = -9.29 , (2) 

rriu = ma 5.5 MeV , m^ = 135.7 MeV , 

which fits TO^, mx, m^i', fm and empirical Alud- 

The thermal part has a coupling to the Polyakov loop 
(spatially homogeneous A4 background) in a form as 



thermal 



-2T 



^ (2^)3 I 



In det 



(2^) 

+ lndct[l + Lte-(^.(f)+A')/7']|, 



l + Le-(^'(f)-^)/^ 
(3) 



where the Polyakov loop is a 3 x 3 matrix in the funda- 
mental representation in color space defined by 



L{x) —V exp 



/3 

-ig I AxiAi{xTX4) 
'0 



(4) 



and in this paper the Polyakov loop sometimes refers to 
the traced one after average, that is. 



1 



(trL) , 



e 



1 



(trLt) 



(5) 



as long as no confusion may arise. In a simple mean-field 
approximation for the Polyakov loop (group) integration, 
the determinant explicitly reads; 



det 



det 



1 + Le-(^-'^)/^ 



1 



-3(6-^)/^ 



+ 3€e-(^-'''/^ + 3^e-2(--'^)/^ 



1 + Ue 



-{E+ll)/T 



1 



-3(e+M)/T 



+ 3^-e-(^+M)/T^3^^-2(e+p)/T_ 



(6) 
(7) 



It is important to note that a positive /i induces £ > £, 
while £ = £ at ^ = 0. The finite-temperature field theory 
tells us that the traced Polyakov loop gives the exponen- 
tial of the free energy cost by a test quark, i.e. £ — gT^'^I^ 
up to normalization (or energy offset), and the anti- 
Polyakov loop by a test anti-quark, i.e £ = e"-''*/-^ [28l |. 



Therefore, /„ — > 00 and thus 



signifies quark con- 



finement, so that the Polyakov loop serves as an order 
parameter for the color deconfinement phase transition. 
In the presence of dynamical quarks in the color funda- 
mental representation, however, neither £ nor £ can be 
strictly zero due to screening. If the color screening is 
stronger, the free energy cost is smaller, and the Polyakov 
loop is larger accordingly. In a medium with > the 
test anti-quark is screened more efficiently than the test 
quark, that means £ > £ [1 [13, IH, ill . 

We shall choose the Polyakov loop potential as [2l| 



Polyakov = -6 •T{54e-'^/^«" 

-|-ln[l - 6^1-3(^^)2 + 4(^3 I 



(8) 



Here, there are two parameters a and h in the above 
ansatz. We fix a = 664 MeV to reproduce ~ 270 MeV 
in the pure gluonic sector and h — 0.03 to yield Tc ~ 
200 MeV for the simultaneous crossovers of deconfine- 
ment and chiral restoration. 



III. RESULTS 

We are now ready to proceed to the numerical cal- 
culations using the PNJL model. Here let us focus 
on the entropy per baryon density ratio, s/ub- This 
is because the adiabatic hydrodynamic expansion con- 
serves s/ub along the time evolution. It is easy to con- 
firm that s/ub is a constant indeed from conservation 
of the entropy current and the baryon current, that is, 
(d/dr)(s/nB) = u^df^{s/nB) ~ readily follows from 
dfi{su^) ~ and d^{nBU^^) — 0. Hence, one does not 
have to integrate the hydrodynamic equation to draw 
the time-evolution path, which is simply inferred from 
an s/riB-constant line . This is the case insofar as the 
expansion is fast enough to make the system thermally 
isolated and the entropy production due to dissipation is 
negligible. 
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FIG. 1: Chiral susceptibility with respect to light (u and d) 
quarks defined by — d'^S7/dm^ in the unit of A as a function 
of and T. 

A. Case without constraint on Us 

We first take a quick look at the results without con- 
straint on Us ■ The quark density is specified by the quark 
chemical potential (or one third of the baryon chemical 
potential) which is common to all three flavors. In the 
results presented in this subsection, thus, the net strange 
quark density rig is nonzero. 

We solve the following gap equations self-consistently; 

dfi _ on dn _dn _ 

d{uu) ~ d{ss) ~"dl ~~d£~ ^' 

with assuming isospin symmetry (dd) = (uu). In this 
way we have the chiral condensates, the Polyakov loop, 
and the anti-Polyakov loop as functions of T and /i. To 
illustrate the phase structure, we show the light-quark 
chiral susceptibility in a density plot in Fig. [1] 

We can perceive from Fig. [T] that the critical region 
extends in the vicinity of the second-order critical point 
located at (T, = (315 MeV, 100 MeV) on top of the 
enhanced strip along the chiral crossover. Because our 
approximation neglects the soft-mode fluctuations, the 
EoS obtained in this work may miss the singular con- 
tribution to thermodynamic quantities near the critical 
point [1^, [lO, [3l[. We already know, however, that 
the PNJL model can reproduce the pressure behavior 
at /X = 0, which implies that the soft-mode contribution 
is not significant there. We can then anticipate that the 
singular contribution would become important only in 
a narrow region surrounding to the critical point (indi- 
cated by the bright colors on the density plot in Fig. [T]) . 
In what follows we shall draw the isentropic trajectories 
onto this phase structure not taking account of soft-mode 
fluctuations. 

By substituting the solution of the gap equations into 
ri, the thermodynamic functions such as the entropy 
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FIG. 2: Isentropic trajectories on the ^-T plane in the case 
without constraint on ris. The (blue) triangle marks the lo- 
cation of the critical point. 



and the baryon number of our current interest obtain 
by means of the thermodynamic relations; 

dn idVL 

'^-df' "B = -3^. (10) 

Figure [2] shows the numerical results for the isentropic 
trajectories for various values of s/ng without imposing 
constraint on n^. 

We note that Fig.[2]is reasonably consistent with Fig. 9 
in Ref. ^| which employs the two-flavor PNJL model. 
Moreover, the trajectories look qualitatively similar to 
the results in the NJL model as in Ref. [2^], while the 
value of s/ns associated with each trajectory is greater 
in our case than in the NJL model study. This has an 
intuitive interpretation. The PNJL model is composed of 
quasi-quarks and a part of gluons, so that s has steeper 
behavior near Tc and grows larger above Tc as compared 
to the NJL model. As for ns, because of the Polyakov 
loop average, ub has steeper behavior as well, but it does 
not exceed the NJL model value. The ratio s/ub in the 
PNJL model, therefore, results in mild sensitivity to the 
steepness of s and ub near Tc, leading to the similar 
trajectory curves to the NJL model results. In contrast, 
the magnitude of s/ub corresponding to the trajectory 
becomes larger in the PNJL model as a consequence of 
the giuon degrees of freedom contributing to s. 

Comparing our results to the lattice QCD simulation 
with the Taylor expansion we see that our estimate 
of s /ub (for instance comparing the s /ub = 40 curve in 
Fig. [2] and the s/n^ = 45 curve in Ref. [l^) improves an 
agreement. 

We remark that all the trajectories in the low density 
side must go to ^ > M„ as T decreases because rts ^ 
when fi < Mu and T ^ 0. In the hadron phase at small 
T and moderate /z, in any case, the PNJL model de- 
scription is not realistic as it cannot describe the nucleon 
and nuclear matter. We should be aware that the PNJL 
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FIG. 3: Isentropic trajectories with Us — imposed, which 
is relevant to the relativistic heavy-ion coUisions. The (blue) 
triangle marks the location of the critical point in this case. 
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FIG. 4: Pressure along the respective isentropic trajectories 
given in Fig. [S] 



model works well near and above Tc and fic but not far 
below them. 



B. Case with constraint on Us 



We could have placed a plot here for the induced /Zs 
as a function of T. For later convenience, however, we 
postpone showing it and let us turn to the validity of the 
Taylor expansion as utilized in the lattice QCD simula- 
tion. 



We next proceed to the case with imposing Ug = to 
emulate the situation in the high-energy heavy-ion colli- 
sions. We should determine self-consistently solving 



together with other gap equations in Eq. 

The phase structure is only slightly changed by the 
constraint. The critical point moves from {jJi^T) = 
(315 MeV, 100 MeV) to (^,T) = (317 MeV, 100 MeV). 
The resultant isentropic trajectory as shown in Fig. [3] 
thus takes a very similar shape to the case without con- 
straint Us = 0. We further calculate the pressure along 
the isentropic trajectories on Fig. [3] and make a plot of 
Fig. 3) Here we did not normalize the pressure by 
as is often the case, for the density contribution cx 
becomes predominant as s/ub goes small. 

From Fig. [3] we see that the strangeness neutrality con- 
dition has a noticeable effect on the trajectory only at 
high T . This is because the s-quark constituent mass, 
Ms, is still heavy around the chiral crossover with re- 
spect to u and d quarks, and so ~ even without 
constraint as long as T is low and /i is smaller than M^. 
If ^ surpasses , we would recognize a difference by the 
neutrality effect in the low-T region, but then, we have 
to consider color superconductivity at such high density, 
which is not within our current scope. 

The isentropic thermodynamics hardly changes, as ob- 
served in Fig. m until s/ns 5~, 20. This behavior agrees 
with the lattice results [l^, [Tsl in which no significant 
s/ub dependence has been found. 



IV. VALIDITY OF THE TAYLOR EXPANSION 

All of our calculations result from the PNJL model 
in the mean-field approximation and do not rely on the 
Taylor expansion unlike the lattice QCD simulation. Al- 
though we have no need to carry the expansion out, it 
is interesting to compare "full" results of our numerical 
calculations and "truncated" ones to verify how nicely 
the Taylor expansion works. 

Within the framework of the PNJL model we have 
solved s and ub as functions of u . Now we shall expand 
them into the series like Ref. [33| as 

n— 

, . ,. ^ d^"+inB(0) 2„+i 

n=0 ^ 

which is to be validated if no singularity associated with 
the first-order phase transition lies along the /x-direction. 
It should be mentioned that the derivative in Eqs. ([72]) 
and (fT5|) is the total derivative in a sense that it acts on 
the implicit /^-dependence in the mean-fields. In this way 
the mixing effect can be correctly taken into account in 
the model treatment j33j . From symmetry s is an even 
function of /i and ub is an odd function. 

Let us explain the numerical procedure in details to 
make a comparison between the results with and without 
truncation of the degree in the Taylor expansion. We first 
approximate s{p) and nsifJ-) by the Taylor expansion 
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FIG. 5: Comparison between the full results (solid curve) and 
the truncated Taylor expansion up to the ^''-term in s{fj,) and 
the n^-term in nsi^^j)- The (green) dotted and (blue) dashed 
curves represent the estimates by the Taylor series fitted to 
reproduce the full results in the ranges from zero to 50 MeV 
and 100 MeV, respectively. Four curves are for s/ub = 40, 
30, 20, 10 from the left to the right. 
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FIG. 6: Comparison between the full results (solid curve) and 
the truncated Taylor expansion up to the ^^-tevra of s{^) and 
the ^^-term of n(/x). The (green) dotted and (blue) dashed 
curves represent the estimates by the Taylor series fitted to 
reproduce the full results in the ranges from zero to 50 MeV 
and 100 MeV, respectively. Four curves are for s/ub = 40, 
30, 20, 10 from the left to the right. 



with sufficiently large number of M and N so that the 
coefficients in the first several terms barely change with 
an increment of M and N . 

We should further specify the fit range of /i e [0, ^o] 
to read the Taylor expansion coefficients. In principle, if 
the "exact" calculation were possible, /xq could be zero 
or there is no /io dependence at all. Even in the mean- 
field level, however, we are far from the exact calculation. 
Here we have chosen /xq = 50 MeV and = 100 MeV. 
In fact, ^0 is a parameter which controls the precision 
in the determination of the Taylor expansion coefficients 
in the same sense as using the multiple-point formula for 
the numerical differentiation. One might have thought 
that /io — 500 MeV, for instance, can cover the whole 
density region in Fig. [31 but such a choice would bring 
artifact from outside the radius of convergence. 

Once we fix the Taylor expansion coefficients of s(/i) 
and nB(/i), then we cut the series at smaller M and N . 
We will elucidate the leading-order case (M — 2 and 
iV = 1) and the next to leading-order case [M — 3 and 
iV = 2) below. 

Figure [5] shows the s/ns-constant curves with and 
without truncation on the ji-T plane. The solid curves 
represent the full results without truncation. We have 
drawn the (green) dotted curves by choosing N = 2 
for s{fi) and M = 1 for nsi^) and the fit range as 
/iQ = 50 MeV. We see that the truncated series can 
approximate the curve for s/ns =40 well, and the curve 
for s/ub — 30 is still close to the full estimate, but the 
dotted curve for s/hb = 20 has a huge deviation from 
the corresponding solid curve. If we determine the Taylor 
expansion coefficients in a wider range as /uo = 100 MeV, 
as shown by the (blue) dashed curves, the agreement be- 
comes better, of course. Not only the curve for s/ns = 40 



but also for s/ns = 30 agrees quite well with the full re- 
sults. Besides, the curve for s/ub = 20 is significantly 
improved at T > 200 MeV, while it does not fit the full 
results at lower T. It seems that the curve for s/ub — 10 
is too far from /io to perceive the effect of changing /io- 

It is intriguing to increase the truncation degree of the 
Taylor expansion to discuss how much the approximation 
is improved. We leave the terms up to the /x^ order in the 
expansion of s(/i) (i.e. M = 3) and the /x^ order in the ex- 
pansion of nB(/i) (i.e. = 2) and have found the curves 
in Fig. [6l Although the results at low T become slightly 
better than Fig. O it is unexpected that the agreement 
at high T goes worse! This undesirable poor convergence 
turns out to stem solely from the Taylor expansion of 
ns(/i)- 

To examine the problem concretely, we present a plot 
for nsifJ.) in Fig. [T] The truncated results at = 1 
(green dotted curve) can well describe the bold solid 
curve which is the full data. As we go to the higher order, 
however, the truncation leads to a larger deviation from 
the full results, implying that the expansion seems to fail. 
If the numerical accuracy is arbitrarily good as needed, or 
equivalently, /xq — > oo, the Taylor expansion would work 
with very small coefficients for the higher order terms in 
the expansion series. In practice, however, the available 
accuracy is limited, and the artifact may enter, which 
eventually goes wrong for /i away from the fitted region. 
Here, at the same time, we should emphasize that there 
is no such weird behavior observed in the expansion of 
s(/x); the higher order we take account of, the better con- 
vergence of s(/x) we can reach, as is naturally expected. 
In view of this, hence, it could be conceivable that there 
may be some profound reason why only the expansion of 
nB(/x) is dangerous. 
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FIG. 7: The bold solid curve represents the full data of nsiiJ.) 
in the model calculation which is Taylor expanded in the fi- 
range of [0, 100 MeV] . The thin solid curves represent the 
truncated results up to A'^ = 1, Af = 2, and = 3 as indi- 
cated by the labels. The dotted curves are the Pade-improved 
resuhs for iV = 2 and = 3. 

We can thus learn the following important lessons from 
these analyses using the model. 

1) The Taylor expansion method does not work in the 
low-T and high-/x region. The adiabatic path for s/n^ < 
20 deduced from the expansion may be totally different 
from the true results. This value of s/hb happens to 
be close to the threshold below which the pressure starts 
moving apart from that at /i = as seen in Fig. [H This 
fact implies that the Taylor expansion may well work 
only in the regime where the density effect is not such 
appreciable. 

2) Even at high T where the Taylor expansion is be- 
lieved to work nicely, the expansion of the baryon density 
ns may be problematic suffering from the uncertainty in 
the higher order terms than 0{fi^). One of the simplest 
remedies for this pathological expansion is the Fade im- 
provement as addressed in Ref. [ll|. We can see in Fig.[7| 
that the Fade-improved results, i?2[l,l] for A^ = 2 and 
i?3[2, 1] for A^ = 3, certainly come close to the full curve. 
[-R2[l,l] denotes ci/x(l -I- C2/Lt^)/(1 -I- ca/z^) with ci, C2, C3 
fixed to yield the original series up to A^ = 2 and i?3[2, 1] 
should be understood likewise.] We note that the entropy 
density s does not have such a kind of expansion problem 
at all. 



V. STRANGE QUARK CHEMICAL 
POTENTIAL 

As promised in Sec. IIIIl the final topic discussed in 
this paper is the induced chemical potential for strange 
quarks to keep the strangeness neutrality. In general the 
condition of ris = requires a positive fXg. This can be 
intuitively understood as follows. 

The strong interaction does not change the number 



FIG. 8: Induced fia necessary to keep Us = for various values 
oi a/uB in the PNJL model calculation. 

of strange quarks but can make strange particles in a 
process, for example, such as 7r~ -I- p ^ K'^ + A and 
7r~ -I- p ^ K'^ -t- Yj~ etc. The strangeness of A and YP , 
is negative one, meaning that they contain positive 
one strange quark. Therefore, if p and n are abundant 
at finite baryon chemical potential, the strong interac- 
tion pushes strange quarks into strange baryons , that 
results in /is > 0. [For more phenomenological details see 
Ref. [Ill ^-iid references therein, and see also Ref. |35].] 

This sort of dynamics is completely missing in the N JL 
model without color confinement. If one solves 71^ = 
in the three-flavor NJL model, it ends up with /^^ = 0, 
which is unphysical. In fact a positive /Xg has been con- 
cluded in Ref. ^Ifli] as it should be. In this section we will 
see that the FNJL model has a crucial advantage in de- 
scribing the induced /is correctly because it encompasses 
the confinement physics. 

Let us first show our numerical results in Fig. [51 A non- 
zero and positive /is certainly appears in the FNJL model 
unlike the NJL model. The numerical values are quali- 
tatively consistent with the lattice results in Ref. [13], 
though the quantitative comparison is not straightfor- 
ward. We would say that this demonstration of /is > 
adds another example to the successful FNJL model ap- 
plications besides the bulk thermodynamics. The rest of 
this section is devoted to explaining how the Polyakov 
loop makes it possible to accommodate the induced /is. 

We recall that the Folyakov loop coupling takes a form 

of 

lndet[l + Le^(^-^)/^] -f Indet [l -I- e-(=+'^)/^] , (14) 

in each flavor sector. As long as fi is small Eq. ([T3]) is 
well approximated as 

c.3e-^/^(£e^/^ + £e-'^/^). (15) 

For light flavors at /i > 0, therefore, the source weight 
for anti-quarks is larger than that for quarks which yields 
£ > ^, as is consistent with the argument given below 



7 



Eq. ([7]). It might be a bit confusing but ie^l'^ is the 
source for anti-quarks because the derivative oi it with 
respect to I gives the anti-Polyakov loop £. Then, in the 
heavy flavor sector, the neutrality condition means in the 
same approximation, 



(TO) log(///) = ^^ 



_^(^e'^=/^ + Ze-^=/n =0, 

Oils 



which immediately leads to 



T. 



= -log(^/^) 



(16) 



(17) 



Interestingly enough, the above equation pT)) holds ap- 
proximately in the entire /i-T plane with only a 3% vio- 
lation at worst! This is an intriguing relation discovered 
in a heuristic manner in the PNJL model and could be 
tested in the future lattice simulation. 

One might wonder why the PNJL model could yield 
a positive /i^ though it does not properly describe the 
confined baryons such as -p and n. As a matter of fact, 
at sufficiently high temperature, the thermal system may 
well consist of mesons and quarks rather than baryons, in 
which TT^ +p JC" + A, for instance, should be replaced 
by TT^ -|- u ^ K'^ -|- s in terms of quarks. The Polyakov 
loop mediates the mesonic correlation through the color 
average, so that this kind of process is to be taken into 
account in the PNJL model. Moreover this process via 
quarks is rather realistic because in the hadronic phase 
the cross section of kaons is small and rig = would no 
longer hold at the later stage of evolution [s^] . 

From Eq. p?)) it is easy to confirm that /is — > when 
T grows large. If T exceeds Tc, both the Polyakov loop 
and the anti-Polyakov loop come close to unity, and thus 
their ratio is nearly one, the logarithm of which is zero. 
This naturally coincides with the intuition that decon- 
fined quarks have no correlation and /is = corresponds 
to ris = just like in the NJL model. In contrast to 
the high-T situation, the confined phase at small T has 
^ ~ ~ 0. The ratio of I and I could be any number. If 
i goes to zero much slower than the ratio can become 
arbitrarily large. 

In the best of our knowledge Eq. p?)) is the very first 
demonstration that the discrepancy between I and i at 
finite /i has a physically significant consequence. The 
importance oi I has been overlooked maybe because 
the difference, £ — £, is negligibly small (only a few % at 
most) as compared to i(f -t-£). The essential point is that 
we consider not the difference but the ratio, which 
may take a huge value when £ ~ and £ ~ 0. 

We plot the induced [is on the [i-T plane in Fig.O The 
functional shape is quite characteristic at small tempera- 
ture. Let us disclose the origin of this triangle-peak struc- 
ture at low T separating the /^--region into three pieces: 

i) Quark-regime — At small \i up to the peak position, 
the induced /is rises linearly along with /i. In this re- 




FIG. 9: Induced \is (or logarithm of the ratio between t and 
€) on the /i-T plane. 



gion the mesonic correlation is the governing dynamics 
as already explained above. 

ii) Diquark-regime — The turning point of /i where 
/is starts decreasing corresponds to the chemical poten- 
tial with which the diquark excitation ie~'^'^'^~>^)/'^ is 
energetically more favorable than the quark excitation 
le~'^'~^^^/'^ . Using Eq. ([T7]) we can derive the condition 
/i, > ^£ for the diquark excitation overcoming the quark 
excitation. Since the constituent quark mass is 336 MeV 
in this model, the threshold should be /i ~ 110 MeV. 
This estimate is really consistent with our numerical re- 
sults shown in Fig. [9l In this //-region /is decreases be- 
cause diquarks behave like anti-quarks in color space. In 
other words the mesonic (quark-anti-quark) correlation 
is taken over by the baryonic (quark-diquark) correla- 
tion which carries nonzero baryon number and so /is is 
partially canceled by this effect. 

iii) Baryonic-regime — If /i is greater than e the color 
singlet contribution, e"^^*^"^-*/^, is dominant, which is 
interpreted as the baryonic excitation. In this regime the 
Polyakov loop is decoupled from the dynamics and there 
is no confinement effect, even though the Polyakov loop 
is zero. (Such a state is recently named the "quarkyonic 
phase" [alll^,!!!!!!].) Then /is = suffices for n, = 0. 



VI. SUMMARY 

We have calculated the isentropic trajectories on the 
phase diagram in the /i-T plane using the PNJL model. 
Our results are quantitatively consistent with the lattice 
results in the high-T and low-/i region where the lattice 
data is available by means of the Taylor expansion. 

To test whether the Taylor expansion is under the- 
oretical control or not, we have expanded our numeri- 
cal results of the entropy s and the baryon density ub 
into polynomial series in terms of /i. We have confirmed 
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that the Taylor expansion can access the trajectories 
for s/ub = 40, 30, and 20, but cannot for s/ns = 10 
at which the isentropic thermodynamics differs substan- 
tially from that at zero density. We have also realized 
that n B (m) has pathological behavior if expanded at high 
temperature, which can be cure by the Pade approxima- 
tion. 

Finally we have discussed the induced strange quark 
chemical potential ^.s to keep rig = in the system ruled 
by the strong interaction. We have found the interesting 
relation between /ig and the ratio of the Polyakov loop t 
and the anti-Polyakov loop This prediction could be 
confirmed in the lattice QCD simulation. 

In this work we neglected the effect of the soft-mode 
fluctuations around the critical point. This is one impor- 
tant direction of the future extension. Another interest- 
ing direction is the origin of the poor convergence in the 
expansion of ub (fJ-) ■ This is not fully understood in the 
present work; other thermodynamic quantities like s(/i) 
have a smooth expansion but only n_B(/i) fails. It would 



be also interesting, if the baryon number susceptibility 
does not have such a problem of expansion (we guess 
so), to validate the idea that the QCD critical point is 
to be deduced from the radius of convergence using the 
model. Although these issues are all beyond our current 
scope, we believe that the present research contributes to 
opening these extensions. 
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